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Abstract 

Spin currents in a two dimensional electron gas with Rashba-type spin orbit coupling are derived 
from a spin connection. Using a functional integral method, we recover the result derived by Sinova 
et aZ.[Phys.Rev.Lett.92, 126603(2004)] for a uniform electric field and in the absence of impurities. 
We extend this result to inhomogeneous electric and magnetic fields. We find that non-uniform 
magnetic fields can give rise to spin currents that are independent of the Rashba coupling and 
hence are less susceptible to impurities than in the case of uniform electric fields. 
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Spin-orbit coupling (SOC) of the conduction electrons in two dimensional systems is of 
special importance in magnetic and semiconductor materials. It was recently realized that 
SOC can be used to manipulate the spin of the conduction electrons in semiconductors. An 
interesting outcome of this coupling is that an equilibrium spin current in semiconductors 
becomes possible due to SOC in the presence of an in-plane electric field" This topic is 
currently of considerable interest since the question on the size of this current is still open. 
One of the reasons for this lack of agreement is that there is no rigorous definition for these 
currents because they are not conserved in a medium with a spin orbit coupling 2 * 3 -. 

In this note, we point out a well defined procedure for calculating spin currents through 
first defining a spin connection 4 similar to the procedure of defining charge currents with 
respect to a U(l) connection. We discuss the simple two dimensional system of an electron 
gas with spin-orbit Rashba-type coupling 5 . The Rashba coupling describes well the dynam- 
ics of conduction electrons in semiconductors, e.g., GaAs, which are potentially important 
materials for spintronics. Similarly, the Rashba model is of interest to conduction electrons 
in magnetic thin films and the damping problem of the magnetization in these materials. 
We will treat two cases; first a 2-dimensional electron gas with SOC in the presence of an 
electric field directed in-plane, and second with an inhomogeneous magnetic field and elec- 
tric field. The first case has already been treated by Sinova et a/.— and others, while the 
non-uniform magnetic field case is new. It also has a universal behavior and hence less 
susceptible to be destroyed by impurity scaatering. 

The recent result by Sinova et al£ that spin currents are possible in systems described 
by Rashba coupling attracted considerable attentio n 7 ^ 9 ! 10 ' 11 ! 12 ' 13 ! 14 ' 15 because of its universal 
features . The value of |e|/87r for the spin Hall conductivity is independent of the strength of 
the SOC. Since spin is not conserved in this system, much debate has been centered around 
the question of what is the correct way to define the corresponding spin current. In the 
following, we argue that the best way to approach the question of what is the appropriate 
equation for the spin current, is by first identifying the associated SU(2) spin connection to 
the current, if such one exists. Then, we calculate the effective action of the 2DEG from 
which a definition of the spin current follows unambiguously by differentiation with respect 
to the SU(2) connection. This procedure is well known in Gauge theor y 16 ' 17 . 
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The Rashba Hamiltonian for a two-dimensional electron gas (2DEG) is given by«^ 

where A is a coupling constant, m is the effective mass of the electron in the lattice, a is a 
Pauli matrix vector, and we use units such that U = c — 1. The spin vector of the electron 
is S = !<7 and its magnetic moment is fi = — ^- (e < 0). The action of this system is 
clearly U(l) invariant and this gives a continuity equation for the charge. Therefore, we 
naturally seek a continuity equation for the spin by first finding under what conditions this 
system has a SU(2) symmetry which is generated by a spin charge. This question was 
first asked by Frohlich and Studer 4 for a general non-relativistic particle interacting in an 
electromagnetic field. They were able to show that such systems exhibit a U(l) x SU(2) 
symmetry and hence a corresponding charge current and spin current follow directly from 
this gauge symmetry. Following this procedure, we can ask similar questions for the Rashba 
system and this will allow us to identify a spin connection and the corresponding covariant 
derivatives. It is well known from Gauge theories, that currents with sources will not be 
continuous but only covariantly continuous. 

In mathematical terms, we want to see if it is possible to express the equation of motion 
for the electrons in the following form 

Bib 1 

{ ~dt = ( VkV ) ^ + eA °^ + 9fiY ° ' SV; ' (2) 

where g = 2 for an electron in vacuum and the spatial covariant derivative is T>k = + 
ieAk + i(g — 1)//Vfc ■ S. A is the usual U(l) field and V is the sought after tensor field 
associated with the SU(2) gauge or spin charge. This way of writing the covariant derivative 
clearly displays that charge currents are associated with the charge e of the electrons while 
spin currents are associated with the magnetic moment \i. Therefore a spin current j k = j^S a 
(summation implied) may exist in a system with zero charge but nonzero magnetic moments 
and will satisfy a 'continuity' equation, T>i~j =0; k = 0,1, 2, 3. This is the result we show 
in this paper within a linear response approach. Hence in the following we neglect second 
order terms in the potential. This is equivalent to taking e and fi as small parameters. 

The case of a constant electric field in the xy-plane and in the x-direction has been well 
studied in the recent literature. Therefore in the following, we focus on the homogeneous 
time-dependent case. The requirement of gauge invariance under U(l) x SU(2) for the 
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Rashba action requires that V = — B and Vk a = — £kap Ep for k = x, y and a, (3 = 1,2,3 
(We use Greek or numerical indexes to denote spin components) . The Rashba coupling A can 
be thought of as due to a fictitious electric field perpendicular to the xy— plane of the 2DEG, 
A = 4^2-^2- The Rashba interaction is therefore replaced by H S oc — ^ (Py&i ~ Px^) — 
-^•p x aj,E x . For an electric field E = —^dA/dt in the x-direction, the interaction term 
becomes to first order in the potential 

g 

H int = — A x (r, t)(p x - Xma 2 ). (3) 
m 

In the following we take E z » E x and hence the additional spin orbit coupling due to 
the E field in the x-direction can be neglected. The E z field will therefore be treated as a 
background field in the following. 

The generating functional of the theory is given by 

Z [A, V] = 

^(t,r)4>(t,r)Ktf(t,r)il>(t,r)]}, (4) 

where we have included scattering due to impurities by a potential U , with (U(r)U(r')) = 
K5 (r — r') where K = 1/ (mr) and r is the scattering. The propagator of the theory is given 
by 

G- 1 = id t - ^vr 2 - [S • V fc , vr fc ] + - g& ■ V„, (5) 

where [, ] , is the anti-commutator operation, 7T; = pi + eA[ and g = 2 for an electron in 
vacuum. In the absence of impurity scattering, the path integral can be easily integrated 
to get the effective action of the theory 

T e// [A, V] = -i Tr logz G~\ (6) 

where the trace is over space and spin indexes. The case of nonzero-potential in the ladder 
approximation has been treated in Ref. 6] 

using a Boltzmann-type approach. To simplify 
our discussion, we set the impurity potential to zero in the rest of this paper since our aim 
here is mainly to demonstrate the spin connection approach. We first briefly treat the case 
of a constant electric field and then discuss the inhomogeneous magnetic and electric fields. 
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The spin current at x = (t, r) in the y-direction and polarized in the ^-direction is therefore 
given byi£ 

1 5T eff 

- TTr 7- \ = -jysi*)- (7) 

This definition reduces to the usual definition of spin currents^. In the linear response 
regime, we need only to calculate G ( x, x') = Go ( x, x') + E ( x, x') to first order in E x 
where Go is the 2x2 Green's function in the absence of the E x field, and its diagonal terms 



are given by G* 1 (w,k) = Gf {u , k) = (u - / (a; - £ ± ie) - A 2 A; 2 , while the off- 
diagonal ones are Gq 2 (u , k) = Gq 1 * (a> , k) = (— iX - ifcj,)) / (uj — || ± ze^ - A 2 /c 2 
The poles of the Green function G give the energy of the spin up and spin down states of 
the system. Since G^ = G^ 2 , only the E-term contributes to the (y3)-component of the 
spin current which becomes 

j y3 = -^ Tr {^P y G° ( Px + \m a 2 ) K G } (8) 

with K = {A x ,p x }. To find the spin current, we go to Euclidean time = it. The 
Fourier transform of the spin current for uniform fields is 

where k± = kp ± mX are the momenta of the two bands and kp is the Fermi momentum, 
Fig. Only states below the Fermi energy are occupied. The static limit is therefore given 
by 



(s-I)tH:, (10) 



E(0) y '%nh 

where we have restored physical units. This result is equal to the one derived by Sinova et 
al& and others^ in the absence of any scattering. In the presence of impurities, it was shown 
111 Refs. Q and Q that no matter how small the scattering by the impurities it cancels the 
spin Hall current. Hence it is natural to ask what happens if the electric or magnetic fields 
are not uniform as is usually the case in real devices. 

In order to shed some light on this question, we extend our treatment to non-uniform 
vector potential fields. We only calculate spin currents in the (y3) direction. We allow for 
both an electric field and a magnetic field to be applied to the Rashba electron gas. The 
vector potential A = (A x , A y , A z ) is allowed to depend on both plane coordinates, x and y. 
To calculate the spin current in the static limit, we use a gradient expansion to the effective 




FIG. 1: Fermi surfaces of the two bands in the Rashba model, k± = kp ± m\. The tangential 
arrows are the spin projections on the k x k y plane. The point k = is a singular point. 

action. In QED 2 +\ such an expansion in momentum space gives rise to the well-known 
Chern-Simons term.™ In our non-relativistic case, a real-space expansion is much easier to 
carry out than in momentum space and this is what we do here. In the following, we keep 
only terms up to second order derivatives in the vector potential. The calculation is carried 
out in the transverse gauge, k • A(k) = 0. The spin current in the (y3) direction has two 
contributions, one coming from the kinetic term and the other from the Zeeman term. The 
kinetic term gives a contribution of the form 

d 2 A x (r,u) 



J a (r, u) = (g - 1) ai (u) A x (r, u) + h (u) 



dx 2 



dy 2 

where the frequency dependent coefficients are given by 

ieX f d 2 ki du 



dydx 



oi (u) = — 



hi (u) = — 



2m J (2tt) 2 2tt 
d 2 k\ du 



-tr {a 3 G° (k 1; u t ) a 2 G° (k x , u t - u) } kf, 



Am J (2tt) 2 2tt 
ieX f d 2 ki du 



Ur {a 3 G° (k ls ut) d 2 G° (k l9 u x - u) } kfkf 



Ur {a 3 G° (ki,^) a 2 d 2 G° (k 1; cu, - cu)} kf, 



Am J (2vr) 2 2tt 



(11) 

(12) 
(13) 



(14) 



ClM = £/ ^^{a^^^d^G^^-u^kf (15) 

The first term, ai(a;), is the original electric field related term found above for the uniform 
electric field and is given (in the imaginary-time) by 

a l(w ) = _^(4mA' (16) 



327rmA 2 
+ojtan~ 1 



2X(k F -m\) 



a; tan 1 



2X(k F + m\) 



Hence in the static limit, u — * 0, this term vanishes and hence we choose not to include the 
electric field contribution below which is the same as in the uniform case. In the non-static 
case time dependent derivatives terms of the electric field also appears in the time-dependent 
spin current. The remaining non-zero cu-dependent terms are too long to be given here. For 
the (y3) direction, the Zeeman contribution is nonzero only for a magnetic field with a 
normal component to the plane, 

J Zeeman (r, u) = ~g^2 ( 17 ) 

xd y B z (r,u;) J |^^Tr {a 3 G° (kx, m) a 3 d ky G° (kx, a* - a,)} k\. 

Summing all these terms and making use of the transverse condition for the gauge field A, 
we get in the static limit the following simple result 

f 2 (r, 0) = d y B z (r, co = 0) . (18) 

This is the main result of this communication which is really the Chern-Simons equivalent 
of the Rashba model and is gauge-invariant. The —1 in (2g — 1) can be traced back to 
the Thomas correction. The y-partial derivative can be understood by going to the rest 
frame of the spin, i.e., zero torque frame, where the SO is seen to give rise to a force in the 
y-direction which is the source of the spin current in a direction perpendicular to the charge 
current. We observe that in non-uniform magnetic fields, the conductivity is proportional to 
the magnetic moment of the conduction electrons and is independent of the Rashba coupling 
as was the case for uniform electric fields. In this reformulation of the Rashba model this 
fact is easily understood in terms of SU(2) gauge invariance. Using Maxwell's equations, it 
can be shown that the spin conductivity is approximately equal to the charge conductivity 



multiplied by the magnetic moment. Moreover, it should be stressed that the constant that 
appear in front of the magnetic field gradient is really a property of the topology of the 
two-dimensional A;— space. 

The fact that a non-uniform magnetic field gives rise to a spin current is not surprising. 
Loss and Goldbart 2 ! were able to show that in a ring with a non-uniform Zeeman term can 
give rise to an equilibrium , or persistent, charge and spin currents. In their system, however, 
the persistent current is due to the phase coherence of the wavefunction. In the presence 
of a non-uniform magnetic field, a Berry's phase is added to the usual Aharonov-Bohm 
phase. Enforcing single- valuedness of the wavefunction results in persistent charge and spin 
currents. In contrast, there is here no real-space topological constraint as is the case for 
a mesoscopic ring. Similarly in Ref. |2l| . it was shown that a persistent spin current is 
possible under the action of an s-d exchange term induced by a non-uniform magnetization. 
The magnetization in that case plays similar role to that of the non-uniform 5-field in the 
Zeeman term. In Fig. 1, the k x = k y = point is a degeneracy point and hence will give 
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rise to a point source term for the Berry curvature, Eq.(9) in Ref. 
we are dealing with a ring structure in k space as opposed to real space in the above two 
examples. 

In summary, we have solved for the spin Hall conductivity using a path integral approach. 
The spin current was defined via a spin 577(2) gauge potential. Our main two results are 
first of all that the conductivity in the static limit for a uniform electric field was found to be 
the same as the one calculated by the Kubo formula. Second, we found that the inclusion 
of an inhomogeneous magnetic field in the z— direction alters this result by adding a new 
component to the spin current. This component is proportional to the classical force on 
a magnetic dipole in a non-uniform magnetic field and is still independent of the Rashba 
coupling. Hence impurity scatterings are unlikely to suppress spin currents generated by 
non- uniform fields. 

Note added - Recently, it has come to our attention that the idea put forward here of 
using the U(l) x 577(2) symmetry as a basis for describing spin currents has been also noted 
by Schmeltzer—. The reader is refrerred to his work for a more detailed discussion of the 
spin current equation. 
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